INCOMPRESSIBLE FLOW IN POROUS MEDIA WITH FRACTIONAL 

DIFFUSION 



ANGEL CASTRO, DIEGO CORDOBA, FRANCISCO GANCEDO AND RAFAEL ORIVE 



Abstract. In this paper we study the heat transfer with a general fractional diffusion term of an 
incompressible fluid in a porous medium governed by Darcy's law. We show formation of singular- 
ities with infinite energy and for finite energy we obtain existence and uniqueness results of strong 
solutions for the sub-critical and critical cases. We prove global existence of weak solutions for dif- 
ferent cases. Moreover, we obtain the decay of the solution in LP , for any p > 2, and the asymptotic 
- - - behavior is shown. Finally, we prove the existence of an attractor in a weak sense and, for the 

' sub-critical dissipative case with a £ (1,2], we obtain the existence of the global attractor for the 

, solutions in the space H" for any s > (N/2) + 1 — a. 

(N 

■ 1. Introduction 

i We use Darcy's law to model the flow velocities, which yields the following relationship between 

the liquid discharge (flux per unit area) v e and the pressure 

Ph ; v^-k{Vp + gjT) , 

I where k is the matrix medium permeabilities in the different directions respectively divided by the 

• viscosity, T is the liquid temperature, g is the acceleration due to gravity and the vector 7 £ M.^ is 

I the last canonical vector cat. While the Navier-Stokes and the Stokes systems are both microscopic 

equations, Darcy's law yields a macroscopic description of a flow in the porous medium [T]. To simplify 
the notation, we consider k = g = 1. 

In this paper we study the transfer of the heat with a general diffusion term in an incompressible 
, flow. The system which we consider is the following (for more details see |18)): 

^ '. 

(1-1) ^+„.VT=-.A"T, 

~^ (1.2) t- = -(Vp + 7T), 

^ . (1-3) diYv - 0, 

I where v > 0, and the operator A" is given by A" = (— A)"/^. We will treat the cases < a < 2 

OP ■ and denote it by DPM. The case a = 1 is called the critical case, the case 1 < a < 2 is sub-critical 

^■"^ ' and the case < a < 1 is super-critical. Roughly speaking, the critical and super-critical cases are 

mathematically harder to deal with than the sub-critical case. 

In [Hj P. Fabric investigates a system of partial differential equations describing the natural convec- 
; ^ ' tion in a porous medium under a gradient of temperature, which is obtained by coupling the energy 

5^ \ equation and the Darcy-Forchheimer equation. He proved existence, uniqueness and regularity of the 

evolution problem as well as the existence of stationary solutions for the two-dimensional case. More- 
over, a regularity theorem is established and a uniform estimate in time of the second-order space 
derivatives of the solutions of the three-dimensional case is given. In [11] the authors consider the 
large-time behavior of solutions to the system 

jvt + V + Vp — Ra*jT = 0, 
divw — 0, 

Tt~ AT + vWT - = 0, 



X 



Date: April 4, 2008. 

Key words and phrases. Flows in porous media 

2000 Mathematics Subject Classification. 76S05, 76B03, 65N06. 

The authors were partially supported by the grant MTM2005-05980 of the MEG (Spain) and S-0505/ESP/0158 of 
the CAM (Spain). The fourth author was partially supported by the grant MTM2005-00714 of the MFC (Spain). 



1 



2 



A. CASTRO, D. CORDOBA, F. GANCEDO AND R. ORIVE 



describing the natural convection in a porous medium filling a bounded domain in M'^. The asymptotic 
behavior of the solutions was studied using the concept of an exponential attractor, i.e. a compact 
finite-dimensional set invariant under the flow associated with the system, and uniformly exponentially 
attracting all the trajectories within a bounded absorbing set. The main results include the existence 
of exponential attractors as well as their strong continuity in a singular (adiabatic) limit 7 — *■ 0. In 
[TTj . using a different method (Galerkin), it was established a global existence and uniqueness result 
for the strong solutions of the three-dimensional Benard convection problem in a porous medium. 
Furthermore, a Gevrey class regularity is obtained for the finite-dimensional attractor of the system. 
Later, in [T^], the authors deduce the x regularity of the attractor. Combining this with a 
Fourier splitting method, they were able to establish the real analyticity of solutions in the attractor. 

More recently the Boussinesq approximation of the equations of coupled heat and fluid flow in a 
porous medium is studied in [5]. This system corresponds to (|l.ip - (|1.3p with a = 2. They showed that 
the corresponding system of partial differential equations possesses a global attractor. They give lower 
and upper bounds of the Hausdorff dimension of the attractor depending on a physical parameter of 
the system, namely the Rayleigh number of the flow. 

Next, we rewrite the system (jl.ip to obtain the velocity in terms of T. The 2D inviscid case is 
shown in [7j. Due to the incompressibility condition, we have that Av — — curl(curl?;). Then by 
computing the curl of the curl of Darcy's law (|1.2p , we get 

d'^T d^T d^T a^T' 



dxidx3 ' 8x28x3 ' dxl 8x2 



Taking the inverse of the Laplacian 

If 1 / 8^T d^T 8^T 8^T 



At: J \x — y\ \dx1dx3' 8x28x3' 8x\ 8x2 



dy 



and integrating by parts we obtain 

|(0,0,T(.,... , 



(1.4) „(a;,i) =-^(0,0,T(a;,i)) + -lpy / K{x - y)T{y,t)dy, xe 



where 



In sections 2 through 4 we consider the case where the spatial domain can be either the whole 
or the torus with periodic boundary condition. 

In section 2 we obtain results of existence of strong solutions of the system (|l.ip - (|1.3I) under the 
hypothesis of regular initial data Tq S H'^ with s > and a € (1, 2]. The case a = 2 was studied also 
in [8]. For the supercritical case a S [0,1), there is global existence for small initial data Tq G i?^ 
with s > N/2 + 1. Also, in the critical case a = 1, the existence of strong solutions is obtained as in 
[31 115j for the critical dissipative quasi-geostrophic equation. 

In section 3 we present results of global existence of weak solutions. We prove a generalization 
of the classical Leray-Prodi-Serrin condition for the uniqueness of the solutions and obtain global 
existence and uniqueness for the subcritical case. In section 4, we obtain the decay of the solutions of 
(IIIlll-IOl) in and for the L^-norms. 

Since we are dealing with a dissipative system, we study in section 5 some attracting properties of 
the solutions of (|l.ip - (|1.3p in with a source term / time independent: 

(1.5) ^-f fVr + i.A"r = /. 



It is easy to see that T, the mean of the solution T of (jl.Sp . satisfies 

d. 



where 



T{t) = j T{x,t)dx, and f = J f{x)dx. 



Therefore, without loss of generality, we can assume that / and T are always mean zero. In particular, 
we prove the existence of a global attractor in the set of the weak solutions with the weak-topology 



POROUS MEDIA WITH FRACTIONAL DIFFUSION 



3 



of L^g^(0, oo, L^(T^)) and a global classic compact attractor, connected and maximal in iJ* with 
s > N/2 + 1 — a in the topology of the strong solutions. 

In section 6 we present results of local existence and blow up of solutions with infinite energy for 
1/ = and for > in the case a — 1,2. 



2. Strong solutions 

Here we show global existence results of the DPM system (|l.ip - ()1.3|) in the sub-critical case. We 
use a maximum principle for the norm of the solutions of DPM, 

(2.1) ||T(t)||Lp < IIToIIlp withl<p<oo, 

which is a consequence of V • w = and the following positivity lemma (see [5D] and [H]): 
Lemma 2.1. For /, A"/ e L'p with < a < 2 and I < p, it is satisfied 

(2.2) J\frhign{f)A"fdx>Q. 

Theorem 2.2. Let Tq G n Lp with s > and N/{a - 1) < p < oo. Then, there exists T G 
C([0,oo);iJ^), solution of DPM with l<a<2. 

Proof. For a solution of DPM we have 

Tt = - div{vT) - vA'^T. 
We use the equality = A{Ri), where Ri is the Riesz transforms (see [131), to get 

i^||A^T|li. A^+^TA^+^-^(i?.(t'.T))dT-i.||A^+tT|ii.. 

Holder inequality and the Calderon-Zygmund inequalities for the Riesz transforms (see [13]) give 

^^IIA^Tlli, < ||A^+tT|U.||A«+i-t(^;T)|U.-j.||A^+tr||i.. 

By the estimate for the operator A** applied to the product of functions (see [H]) for s > 

(2.3) ||A^(/.g)||L^<C(||/||i,n|A'5lU,|| + ||5lL,'||AVIlL0 K r < q' < ^, - = - + ^, 

r q q' 

we have for q' — p, 

Ij^W^^nh < C||A^+tr|U.(||t;|U.||A^-+^-^T|U, + \\T\\l4^^+'-^v\\l,) - iy\\A'+^T\\l., 

with (1/p) + (l/(?) — (1/2). Now, since v satisfies (|1.4p . again we apply the Calderon-Zygmund 
inequalites obtaining 

\\v\\lp < C\\T\\lp, \\A'+^-^v\\l, < C\\A'+'-^T\\l., 

and (|2.ip gives 

^^IIA^Tiii. < c||rol|L.||A^+^r|U.||A^+i-tr|U, - i.||A^+tT||i.. 

The inequality for the Riesz potential (see [53]) 

(2.4) ||/M/)||l, <C||/||„ 0<(3<N, l<r<g<oo, - = --^, = A"", 

q r N 

for r = 2 and f3 — N/p, yields 

\\A'+^-^T\\l. < C\\A'+^-'i+I^T\\L2. 
We take p > N/{a — 1), and therefore 1 -I- /3 < a, so that 

IIA^+i-t+'^TlU^ < ||A^+^T||2.||A^r||^I^ 
with 7 = {2~ a + 2(3) /a < 1. Applying the last two inequalities we obtain 
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:\\K^T\\l. <C{u^\T,U,)WT\\l.. 



and Young's inequality gives 

ld_ 
2di' 

Furthermore, we have 

||A^T|U.(t) < WA^Toh^e^'K 
From this a priori inequahty together with the energy estimates argument we can conclude the global 
existence result. ■ 

Theorem 2.3. Let < a < 1 be given and assume that Tq £ H'^ , s > {N — a)/2 + 1. Then there is 
a time r = r(||A*ro||) so that there exists a unique solution to DPM with T G C{[0,t), H'^). 

Proof. Since the fluid is incompressible we have for s > (N — a)/2 + 1 
i AllATIIIs = - y" A'TA'{vVT) dx - u\\A'+'^T\\l2 

= - j A"T(A"(wVr) - ?;A'*(Vr)) dx - A'^+trji^ , 

< C||A'^T||L2||A''-(t'VT) - vA''{VT)\\l2 - vWA'^+'^TWl^. 
Using the following estimate (see |13j ) 

||A^(/5) - /A^(.9)|Up < C (||V/|U^|1A«-i.9|Up + |lAV||Lp||g||L~) K p < oo, 
we obtain for p — 2 

ld_ 
2^ 

Applying Sobolev estimates we get 

^^IIA^Tlli. < C(||To|U. + ||A^/2+i+-T|U.)l|A^T||i. - >y\\A^+^T\\l, 
and taking e = s+ ^ — ^ — lit follows 

(2.5) ^I'l^'^ll'^^ ^ + l)(rolli. + |lA^r|liOI|A^T|li.. 

Local existence is a consequence of the above a priori inequality. 

Let us consider two solution and of DPM with velocity and respectively, and equal to 
the initial datum T^{x, 0) = T^{x, 0) = Toix). If we denote T = - T'^ a.nd v = - v^, we have 

ljt\\nh<- jrvVT^dx-vWA^ni.. 

For a = Calderon Zygmund and Sobolev estimates give 

~\\nl^ < C\\T\\U\VT^U^. < C\\T\\U\T'\\h^. 
Inequality (|2.5p implies that ||r^||/f3(t) is locally bounded. Furthermore, we can conclude 

mUt) < ||T||i.(0)cxp(C f \\T'\\H^{cT)da), 

Jo 

which yields uniqueness. 

The case a > is treated differently, we have 

l^^Wnh < ||T|U.niL,||vri|U. - HIAtTiii., 

with q = 2N/{N - a), and p = 2N/a. Since ([^ we obtain 

\\t\\l. < c||Atr|U., ||vt1|Up < cwA^+f-^r^u^ < c\\t'\\h^ 

and finally 

il||r||i.<^i||r||i,||T.||^. . 
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Remark 2.4. For the supercritical cases (0 < a < I), we have the same criterion as 7 for the 
formation of singularities in finite time. In fact, we have that T € C([0, r]; H'^) with s > N/2 + 1 for 
any t > if and only if 



f \\VT\\BM0{t)dt <(X,. 
JQ 



For small initial data we obtain the following global existence result for the supercritical case. 
Theorem 2.5. Let > 0, < a < I, and the initial datum satisfies the smallness assumption 

llrollff. s>N/2 + l, 

for C a fixed constant. Then, there exists a unique solution of (|l.ip - (|1.3p in C([0, oo); i?'*). 
Proof. We multiply by K^^T and, by the Sobolev embedding, we get 

IjtWnl^ < CiWVvU^ + ||VT|U.)||A^r||i. - u\\A^+^T\\h 
< Ci\\T\\L2 + ||A^r|U.)l|A^T||i. - i/||A^+tT||i. 

with s > 2. Thus, we have 



Since 
we obtain 



||A"T||2, < ||Atr|||2 + ||A-'+tr|||2, 



by the assumption of the smallness of the initial datum. ■ 

In the critical case a = 1, we state the following regularity result. 

Theorem 2.6. Let T be a solution to the system ()l.l[) - (ll.3p . Then T verifies the level set energy 
inequalities, i.e., for every A > 

J Tl{t2,x)dx + J J \A^/^Tx\^dxdt< J Tl{ti,x)dx, 0<ti<t2 

where T\ = {T — X)-^-. It yields that for every tg > there exists 7 > such that T is bounded in 
CT([to,oo) X M^). 

With this result we show that the solutions of the diffusive porous medium with initial data 
and critical diffusion (— A)^/^ are locally smooth for any space dimension. The proof is analogous 
to the critical dissipative quasi-geostrophic equation that it is shown in [3]. Analogous result can be 
obtained using the ideas of [T5] to show that the solutions in 2D with periodic C°° data remain C°° 
for all time. 



3. Weak solutions 

In this section we prove the global existence of weak solutions for DPM with < a < 2. First we 
give the definition of weak solution. 

Definition 3.1. The active escalar T{x,t) is a weak solution of DPM if for any (p G Cj?°([0,T] x M) 
with ip{x, t) — 0, it follows: 



(3.1) 0= J To{x)ip{x,0)dx + J j T{x,t){dtip{x,t)+v{x,t)-Vip{x,t)-vA"Lp{x,t))dxdt, 

R" RN 

where the velocity v satisfies (|1.3p and is given by ()1.2p . 

An analogous definition is considered in the periodic setting taking ip £ C°°{[0,t] x T). 
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Theorem 3.2. Suppose Tq £ L^(R^) and < a < 2. Then, for any t > 0, there exists at least one 
weak solution T E C{[Q,t]; L^{R^)) n L^{[0,t]; H'^/'^{R^)) to the DPM equation. 

Proof. To prove the theorem we modify the system (|l.ip - (|1.3p with a small viscosity term and we 
regularize the initial data. In particular, for e > 0, we consider the family Tg of solutions given by the 
system 

BT 
at 

(3.2) Ve = - (Vpe + iTe) , 

divide — 
T,(x,O) = 0e*To, 

where * denotes the convolution, (j)e{x) — e^^(f>{x/e) and 

0eCf=(R^), 0>O, [(j){x)dx^l. 



As we show in the previous section there is a global solution of (|3.2p with E C{[0, t]; H'^ (R^)) for 
any s > 0. We multiply by to get 

i|(||r,||i.) + H|A"/^r,||i.<o, 

and integrating in time 

(3.3) m{r)\\h + 2^ r \\A"^'T,is)\\hds < \\To\\h Vr. 

Jo 

In particular we find 

(3.4) T,eC([0,r];L2(]R^)) and max ||r,(i)||i. < ||ro||i.. 

0<t<T 

We pass to the limit using the Aubin-Lions compactness lemma (see [1 5j ) : 

Lemma 3.3. Let {f^{t)} be a sequence in C{[0,t]; H'' (R^)) such that 
i) max{||/e(i)lk= : 0<t<T}<C 

a) for any Lp E C^(R^), \}pfe\ is uniformly Lipschitz in the interval of time [0, r] with respect 
to the space H'^{R^) with r < s, i.e., 

\\ipfe{t2) - Vfe{ti)\\H'' < Cs\t2 - ti\ < ti,t2 < T. 

Then, there exists a subsequence {fej{t)} and f E C{[Q,t]; H''{R^)) such that for all X E {r,s) 

max \\(pfej{t) - ipf{t)\\H^ ^0 as j ^ oo. 

First, by ([331) we get E C{[0,t]; L^{R^)) and i) in the space i^(R^). Next, we prove that the 
family T^ is Lipschitz in some space H^'''{R^) with r > A^/2 + 2. Since is a strong solution of 
and continuous it follows 

*^ d 
ip—T^dt 



dt 



< max {A{t)}\t2 - til 



fj-^ tl<t<*2 



(3.5) yT,it2)-ipniti)\\H-r -~ 

where 

A{t) - II div(^,;,T,)||^-. + II div{^)v,T,\\H-^ + ,y\\^A"T,\\H-. + e\\^AT,\\H-^. 
Applying the property that the Fourier transform of the product is the convolution of the respective 
Fourier transforms, we have 



m\^-rmi-v)dv <C / (|^r + |r;r)|^(r;)||Tae-?7)M??<(l + kr)||^lk=||Te|U., 
and it yields 

||^A"T,||h-,- < C(^)||r,|U.( J i^-^l^de) < C{r,^)\\Toh^. 

Analogously, 



|^Ar,||H-<c(r,^)||ro|| 



L2. 
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We have 

||div(^)z;,r,||^-,. <C(r)|| di^(^,r,|U^< C(r)|| dh^)|Ui II^T^Hi^ < C(r, (^)||«,|U2 ||r,|U2, 
and by p.4|) and the fact that the velocity satisfies (|1.4p . it foUows: 

||div(^KT,||H- <C(r,^)||ro||i.. 

In a similar way 

II diy{ipv,T,)\\H-r < \\ipv,T,\\H.-r < C(s,(p)||ro||i.. 

From p.Sp . the condition ii) of the Aubin-Lions lemma is satisfied. Therefore, there exists a 
subsequence and a function T E C([0, r]; L^(R^)) such that 

(3.6) T. T in a.e. t and ma.x \\ipTe{t) - ipT{t)\\H>^ ^ asAe(-r,0). 

0<t<T 

We pass to the limit in the weak formulation of the problem (|3.2p . i.e., 



0= J T^{x,0)ip{x,0)dx + J J T,{x,t) {dtip{x,t) + v^{x,t) ■\Iif{x,t) - vK"if{x,t) + el^ip{x,t)) dxdt, 
and we obtain 

T T 

= y To(a;)(p(a;, 0)da; + j j T{x, t) {dtip{x, t) - uK°'ip{x, t)) dxdt + 1™ y j T,{v, ■ ^Lp)dxdt. 
Next, we decompose the non- linear term 

T T T 

J J T,{v, ■ \7ip)dxdt = j j {Te- T){v, ■ '^(p)dxdt + J J T{ve ■ \/f )dxdt. 

K« R« M« 

Using the Fourier transform in the first term we have 
(Te - T){v' ■ W(p)dxdt 



< J \\Ve\\H^Mm-T)Vip\\H-.,2dt 


T 

< max ||(r,~T)V^||^_/. / {\\v,U^ + WA^/'v^U^) dt. 



Due to dnSl) and ((231) we get 

r 

{\\v,\\L2 + \\A'-/^V,\\L2)dt<c{T)\\To\\L2. 



Then, by p.6p we have 

T r 

lim J J T^{v^ ■ 'Vip)dxdt ~ J J -^(^ ' '^f)dxdt 

R« R« 

and we conclude the proof of theorem [ 



Remark 3.4. Analogous result of theorem \3.'2\ follows, with a similar argument in the torus , with 
periodic boundary conditions. 

We continue this section mentioning the existence result of weak solutions obtained for the non- 
homogeneous equation p.Sp . 

Theorem 3.5. Let t > be arbitrary. Then for every Tq E and f E L^{0;T; H^^), there exists 
a weak solution of IfT^) satisfying T{x,0) = To{x) and T E C([0,t];L2) n L'^{[0^t];H"/'^). 
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The proof is similar to theorem 13.21 

Although weak solutions may not be unique, there is at most one solution in the class of "strong" 
solutions in the sub-critical case. This fact is well known for the quasi-geostrophic equation (see (Ij) 
and it is a generalization of the classical Leray-Prodi-Serrin condition, related to the uniqueness of 
the solutions to the 3D Navier-Stokes equation (see [25]). 

Theorem 3.6. Assume that a e (1,2], r > and T a weak solution of DPM with Tq G L^. Then, 
there is an unique weak solution satisfying: 

(3.7) T e C([0, r]; L2) n ^^([o, r]; t) n LP{[0, r]; L«), 

for q > N/{a — 1), and p — a/{a — N/q — 1). 

Proof. We take the difference T ^ — T"^ oi two solutions and T'^ of DPM with same initial 
data. Considering v = — v^, with and being the velocities corresponding to and T^, then 
T satisfies 

BT 

— + • VT^ + • VT + i/A"r = 0, 
at 

or analogously 

BT 

— + dw(vT^)+ div(u2r) + i^A"T = 0. 
at 

We multiply the equation by K^^T and integrate by parts in the nonlinear terms to obtain 



-||A-^T||i. + H|A^(A-^r)||i.< 



{Tv2) ■ {W{A-^T))dx 



(Tiv) ■ {W{A-^T))dx 



We take (l/q) + (2/p) = 1 and it yields 

/ {Tv^) ■ (yiA~^T))dx 

JR 

and 

J (T^) ■ (V(A-ir))dx 
Since V(A"^) = (i?i,i?2) we obtain 



<\\v'\\l4t\\l4^{^-'t)u^ 

< \\T'\\L4v\\Lp\\y{A-'T)\\LP. 



(Tv^) ■ {V{A-^T))dx 



{T\) ■ {■V{A-^T))dx 



<C{\\T4L. + \\T'h.)\\T\\%. 



The inequality for the Riesz potential (|2.4p gives 

\\T\\lp < C\\A^T\\l2 = C\\A^+Ha-^T)\\l.. 
For q large enough, we can get {N/2q) + 1/2 < a/2, and the following interpolation inequality 

||A7IU^<ll/llI.||AVIli;'^ 

with s<r, 0<7<1, yields 



||A#+^(A-^T)|U. < ||A-^T||l,||At(A-^T)||]^I^ 



for J — {a — N/q — l)/a. Finally 
d 



dt 



and 



||A-^T||i.+j.||At(A-^T)||i. <C(||Ti||i, + ||T2||i,)l|A-^r||2l||At(A-^T)||^i^-^) 



-^||A-^r||i. < -dlTilU, + llT^IUO'/^IIA-^Tlli., 
dt V 



which allows to conclude the proof. ■ 

Remark 3.7. // we take Tq ^ L'^ L'' for q > N/ (a — 1) we can construct, as before, a solution that 
satisfies ||r||L9 < ||To||l9, and in particular we have 

Te L°°([0,r];L«). 

Then this solution is unique in this space. 
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Remark 3.8. Suppose that l<a<2,i^>0,s>0, feLPf] i?""^(T^) and 

To e i/"niP(T^), w/iere < i < 

Then, there is a weak solution T of p.Sp such that 

T e C([0, T];H'iT^)) n L2(0, r; i/^+t (T^)). 

T/ie proof follows applying an analogous analysis as in the proof of theorem \2.Si Moreover, as in 
theorem \S.6[ the uniqueness of weak solutions for the non-homogeneous equation (II. 5p is also obtained. 



4. Decay estimates 

Here, we obtain the decay of the solutions of (|l.ip - (|1.3p . The key of the argument is the following 
positivity lemma. 

Lemma 4.1. Suposse a G [0, 2], n = K^, and T, A"T e LP{n) where p>2. Then 

J \T\P-^TA"Tdx>- J (At|r|^)^dx. 
n n 

This lemma is a consequence of different versions of the positivity lemma obtained in [201 13 El [12] • 
The immediate consequence of the previous lemma is the following decay results in the space 
for solutions of (fl~T|) - (fT3|) : 

Corollary 4.2. Suppose that Tq e where p G [2, +oo) and T is a weak solution of JJ.-?|) - (|1.3p . 

(1) If ~ and the mean value of Tq is zero, then 

qe[i,p\, 



\T{t)\\L. < llTollLpexp 

V P 

where Ai > is the first positive eigenvalue of A. 
(2) If n = R'^, then 



iv(|j-2)||roir<''-'' 



lL2 



Nip-2) 
2pa 



with c depending on a and N . 
Proof. We multiply the equation (jl.ip by \T\p~'^T and applying the lemma |4T] we obtain 



dt 



llTWIli, < -2^ / |A^r«pdx. 



In the case $1 = T we get 



d 
dt 



|T(t)||ip + 2i.A?||rW||i, <0, 



and the above inequality gives the exponential decay of ||T(i)||Lg for q G 
In the case il — M.^ , using Gagliardo-Nirenberg inequality, we have 



^||T(i)||i.<-2^c( / \T\^dx 



d 
dt 

with c depending on a and N. By interpolation we get 



T—dx 



with 



Nip-2) 
N{p ~ 2) + 2a 
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Therefore, 

iiT('/^iip__ -J- '?r7/iiTii':;~ ~ 

dt 

since (3 £ (0, 1) and \\T(t)\\L2 < ||To||l2 yields 



T(t)||^, + 2cH|r||^;^||r||f,<o, 



We integrate 



-|||rwr^, + 2.c^%-<o. 



Again, since (3 £ (0, 1), we have 

V roliLP J - P iirjit-^ 



hence 

/3 



" wnwi.' J 

Then, by definition of /?, it foUows the polynomial decay. ■ 

Remark 4.3. As a consequence of the previous lemma for the case ft = M.^ , we obtain the following 
estimate for the L°° -norm: 

\\T{t)\\L^ < \\To\\l^, 

which can be improved as in [5] 

imOlU- < llTolU^ (1 + actWToWl^yi . 



5. Long time behavior 

In this section we study some attracting properties of the solutions of (II. Sp with a E (1,2]. 

We introduce an abstract framework for studying the asymptotic behavior of this system with 
respect to two topologies, weak and strong, depending on the uniqueness of the solution. Each such 
system possesses a global attractor in the weak topology, but not necessarily in the strong topology, 
and in general, they are different. 

First, we recall some definitions of [26j . Let X be a complete metric space. A semiflow on X, 
uj : [0, oo) X X X,is defined to be a mapping ui{t, x) = S{t)x that satisfies the following conditions: 
S{0)x = X for all x £ X; to is continuous; the semigroup condition, i.e., S{s)S(t)x ~ S{s + t)x for all 
s,t > and x £ X, is valid. 

A semifiow is point dissipative if there exists a bounded set B G X such that for any x £ X there 
is a time t{x) such that uj(t,x) £ B for all t > t{x). In this case B is referred as an an absorbing set 
for the semifiow S{t). 

A semifiow is compact if for any bounded B <Z X and t > 0, S{t)B lies in a compact subset in X . 

^ C X is a global attractor if it satisfies the following conditions: A is nonempty, invariant and 
compact; A possesses an open neighborhood U such that, for every initial data uq in U, S{t)uo 
converges to as t ^ oo: 

dist(5(t)'U07-4) as i — ^ cxD. 
Recall that the distance of a point to a set is defined by 

d{x,A) = inf d(a;,?;). 

yeA 

Now, we state the following result about the theory of global attractors (see [26]): 
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Theorem 5.1. Let S(t) be a point dissipative, compact semiflow on a complete metric space X . Then 
S{t) has a global attractor in X. 

Furthermore, the global attractor attracts all bounded sets in X , is the maximal bounded absorbing 
set and minimal invariant set for the inclusion relation. 

Assuming in addition that X is a Banach space, hi is convex and for any x € X , S{t)x : K_|_ —> X 
is continuous. Then, A is also connected. 

In the case that U = X, A is called the global attractor of the semigroup {S(t)}t>o in X. 

5.1. Strong attractor. From Remark 13.81 we see immediately that for any 1 < a < 2 and with 
s > {N/2) + 1 — a, the solution operator of the porous medium equation (|1.5p well defines a semigroup 
in the space iJ". 

We begin this section with some useful a priori estimates of the solutions (|1.5p with f E L^. 

Lemma 5.2. Let T ~ T{x, t) be a solution of lil.5]) . the initial data Tq £ L^ with zero mean value 
and p > 2. Then, ||r||Lp is uniformly bounded with respect to UToIIlp- particular. 



(5.1) 



imh. < {\\to\\l. ^ii/iu^) exp {-^4 + 



and there exists an absorbing ball in L^ . Moreover, for Tq E L^, we get 

(5.2) z^^*""' ||AtT(.)||i.ds < (^llTolli. - Mii^ cxp{-i.A?0 + Mid + i||A-f 

Proof. We multiply the equation ()1.5|) by p\T\P^'^T with p > 2. Integrating, using lemma \4A\ and 
applying Holder's inequality, we get 

(5.3) ^^\\Tit)\\l, + 2u\\A^Ti\\l.<p\\fh4T\\l-\ 
We denote Ai as the first eigenvalue of A. Since T is mean zero, we have 

I1A^t5||^, >A?l|r||^,. 

Therefore, 

|||r(t)||i. + ^||r|U. <||/IU., 



and integrating we prove the estimate ([F 

Now, to prove (|5.2p . we multiply the equation (11.5|1 by 2T and by 

|||r(0||i.+2H|AtT||i. <2||/|U.||T|U.. 
Integrating and applying the Young inequality, we get 

j^Wmwh + HIA^Tiii. < l\\A-^f\\h < 

Integrating and using that Ai is the first eigenvalue of A, we obtain 
(5.4) ||T(i)||i.< ('||ro||i.-MEi')exp{-.A?0+''^'''^ 



On the other hand, integrating between t and t + I, we have 

(5.5) mt + l)||i. + \\A^T{^)\\hds < \\Tit)\\l. + -llA-t/IU., 



and we get the estimate (|5.2p by (|5.4p . ■ 

Lemma 5.3. Let T — T{x,t) be a solution of \1.5\) . Tq G with zero mean value and s > {N/2) + 
1 — a. Then, ||A''T||^2 is uniformly bounded with respect to ||A*To||i2 and there exists an absorbing 
ball in the space H"^ . Moreover, we have 

(5.6) / ||A"+tT||2,dt < +CX) 
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and 

ct+i 



(5.7) V J ||A'' 2y||^2 is uniformly bounded with respect to \\A^To\\l^ ■ 

Proof. We have that a e (1, 2) and that TqG H', where s > {N/2) + l-a. If s e {{N/2) + l-a,N), 
let r = s. If s e [N, +oo), let r be any real number in ((iV/2) + 1 - a, A^). Then, Q W d LP, 

where 

1 1 r a — 1 

2^ N N 

We multiply the equation (|1.5p . with an initial data Tq belonging to LP n by A^"T with 
< a < s. By an analogous analysis as in the proof of theorem 12.21 we get 



^||A'^T(i)||i. + i.||A''+tr||i. < i||A«-t/||i. + c||T||i.||A'^+^+^r||i.. 
Then, by T e L°°(0, oo; LP) and 

Now, using the Gagliardo-Nircnberg and Holder inequalities, we have 



and 



C\\A_^+PT\\l2 < pA''+^T||^, + ^IIA'^Tii^ 



l\\A^T{t)\\l. + ^IIA'^+^Tlli. < i||A'^-f /Hi. + ^IIA'^Tlli. 



Next, following from the above inequality and (|5.2p the uniform boundedness of ||A°T(t)||i2 with 
respect to ||A°To||/,2 can be obtained for a < a/2 from applying the Uniform Gronwall lemma (see 
Remark 15. 4p . 



||A'^T(t + l)||i. < -e^ ((wToWh - ^) exp{-i.A?i} + Mild + i||A-t/||2, + ||A- 

This estimate can assure us that also it gives us an absorbing ball of the solutions in the space TJ" 
with < a < a/2. 

Moreover, we have (|5.6p and (|5.7p for < a < a/2. Therefore, with these estimates and a 
bootstrapping argument, the uniform boundedness of |jAT(i)||^2 is indeed valid for any s > {N/2) + 
1 — a by using the Uniform Gronwall lemma again. This also gives us, as before, an absorbing set in 
the space LL' for any s > {N/2) + 1 — a and the estimates (|5.6p and (15. 7p . ■ 



Remark 5.4. (Uniform Gronwall lemma) Let g, h and y be non-negative locally integrable functions 
on (to,oo) such that 

^ < gy + h, Vt > to, 
at 

and 

g{s)ds<ci, J h{s)ds<C2, J y{s)ds < c^, > to, 
where r, ci , C2 and C3 are positive constants. Then, 

y(t + r) < (^+C2)e°S yt>to. 

The proof of this estimate is shown in [26j . 

Now, we prove a condition to apply the theorem 15. II the continuity of the solutions of (|1.5p in the 
space if* with respect to t. 

Lemma 5.5. Let T — T{x,t) be a solution of Tq e with zero mean value and s > {N/2) + 

I -a. Then, A^T e C{0,t;L^). 
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Proof. By lemma [5731 we have that A'^T £ L^(0, r; _ff 2 ). According to the Aubin-Lions compactness 
results (see [52]) we just need to show that A'^Tt £ L?{Q,t]H~^). From the equation p.Sp we get 

(5.8) l|A^Tt||^_^ < ||Ai+^-t(T«)|U. + ||A^+tT|U. + ||A^-t/||^,. 

Using (|2.3p and the integral formulation of the velocity we have 



|Ai+''-t(Tw)||i2 < C||r||Lp||Ai+'^-tT| 



L9 



where (l/p)+(l/q)=l/2. Now, as in the proof of lemma [5?3l we take r < s such that Tq £ C 
H'' CLP, with 

1 1 r a ~ 1 

2^ N N 

Then, considering q = N/r, T £ L°°{0, 00; L^). Since r > {N/2) + 1 — a (see the proof of lemma 
we have that 

_ N 2N _ , 

r "^iV + 2-2a~^' 

Therefore, 

|jAi+''-fr||i, <c\\A^+'-'^T\\l,' < C||A'^+tT||i2, 
applying the Gagliardo-Nirenberg inequality. Thus, coming to ()5.4p . we get 

||A'T,||^-= < iC\\T\\L, + 1)\\A'+^T\\l2 + \\A'-^f\\L2. 
Finally, by (|5.ip and (|5.6p . we obtain 



\A'Tt{t)\\^-^dt<(x^, 



and we conclude the proof. 



Lemma 5.6. Assume that the initial data, of a solution of equation il.5]) . belongs to iJ* with zero 
mean value and s > N/2 + 1 — a. Then, for any fixed t > 0, the solution operator S{t) is a continuous 
map from H'^ into itself. 

Proof. We consider two solutions T*^^) and T^^^ of the porous medium equation p.5p with two initial 
data T^^^ and T^^^ and velocities w^^^ and v^'^\ respectively. Let T = T^^) - T^^) ^nd v = w^^^ - u^^). 
Then, since div(t)) — 0, we have for any Lp £ that 

(5.9) (rt,(/3) + i/(A^T,A^V') = -{v ■ VT^^\ip) ~ (w^^^ VT, (/?). 

Setting Lp — T , using that div(u(^)) = and Young inequality, we get 

i|||T||i.+H|AtT||i. <||AT(2)|U,.||T||i,, 
such that (l/fZi) + (2/(72) = 1- By Gagliardo-Nirenberg inequality, we obtain 



with a = N/(qia), where we will be choosing qi such that a £ (0, 1). We use again Young inequality 
and we have 

1 d 

2dt"^ ' 2' 

denoting q — 1/(1 — a). Thus, by the Gronwall inequality, it follows 



nii^ + ^IIA^Tlli. <c(i.)||AT(2)||i,,||T||i., 



||r(t)||i. <C(^)||r«-Tf lli.expf / \\AT('\s)\\Uds 



If s e i{N/2) + 1 - a, {N/2) + 1 - (a/2)), then we take r ^ s. If s £ [{N/2) + 1 - (a/2), +00), we take 
r any number in {{N/2) + 1 - a, {N/2) + 1 - (a/2)). Then C iJ^ We choose 

2A^ 

2 + N-2r-a > ^' 

then 

2 + 7V-2r-a 1, 
a = — £ (0, -) and q < 2. 
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Therefore, using the foUowing Sobolev inclusions 

L«(0, r; W^^"') C L«(0, r; H'^+t) c ^^(0, r; c L2(0, r; H'+''- 

we conclude that 

" \\A^T{t)\\Ut < C{T('\r)\\4'^ - Tf Hi. 



Thus, using the Riesz lemma, it is immediate that the solution operator S(t) is a continuous map 
from H'' into itself when s e {{N/2) + 1 - a, a/2]. 

We finish the proof studying the case s > a/2. We do so by checking directly the Lipschitz 
continuity of the solution operator in the space . We consider ip = A^^T in ()5.9p . then 



(5.10) ^I'l^'^ll'^ + ^W^'^^TWl^ = ~{v ■ VT(2), A^^T) - (t;(i)vr, A^^T). 

We estimate the two terms on the right-hand side of the variational formula separately. First, we get 

|(wVr(2), A^^T)! < c{iy)\\A''^ {v ■ VT^^))!]!, + '^\\A'+^T\\l2. 
Using similar estimate as (|2.3p of Kenig, Ponce and Vega (see [Hj), we have 

l|A^-t(«(i)vr)|U2 < ||A^-tT||LPi||Ar(2)||ip, + ||r|iL,i||A^+i-tr(2)||i,, 

with (1/pi) + {I/P2) - 1/2 and (1/qi) + (l/ga) = 1/2. We select 

27V 2N N 2N 
Pi = T7 > -P2 = , qi = 7, 92 = 



N — a a a — 1 N + 2 — 2a 

and, using the Sobolev inequalities yields the following estimate 

(5.11) |(t>(i)vr, a2^t)| < c(i.)||A*r||2,||A^+tr(2)||2, + ^\\a'+^t\\12. 

We estimate the other term on the right side of (|5.10p . Since (w^^^ • VA^T, A*T) = 0, A* and V are 
commutable we have 

|(i;(i)vr,A2^T)| < |lA^(i-(i) • VT)-i;(i) • (A''Vr)||i2|lA''T|li2. 
Using estimate of Kenig, Ponce and Vega (see [I^) we have 

|jA^(i;W . Vr)-i;W •(A^VT)|U2 < \\Av('^\\l.,\\A'T\\l.2 + \\A'v^'^\\l.,\\AT\\l.,, 
with (1/pi) + {I/P2) - 1/2 and (1/gi) + (l/ga) = 1/2. We take 

2N 2N 



a 1\ ~ a 

and using the Sobolev inequalities we get 

|(w(i)vr,A2^r)| < ||A'*+tT(i)||i2||A''+tr|ii2||A'*r||i,2 

< c(z/)||A^+tr(i)||2,||A«T||2, + ^||A^+tT||i2. 
Therefore, considering this estimate and ()5.1ip in ()5.10p . we obtain 

^||A^T||i2 + uW+^T\\l. < c{v) (||A^+f T(i)||i2 + |1A''+Tr(2)||i,) Wnl.. 
So, by Gronwall's lemma and since 

{\\A'+'iT^^\t)\\l2 + ||A''+tr(2)(t)||2,^ dt < 00, 
we get 

||A^r(t)|U2 < C(i., ||A^+^r(|^)|U2, ||A'^+tTo(^)|U2)||A'^(To(^) - Tf )|U2 
and conclude the proof of the lemma. ■ 

Finally, we present the existence of the global classic attractor. 

Theorem 5.7. Let a £ (1,2], > 0, s > (A^/2) + l — a and f E H^^°'r\U' time-independent external 
source. Then, the operator S, such that S{t)To — T{t) for any i > and T solution 0/ (jl.5p . defines 
a semigroup in the space H'^ and satisfies: 

i) For any t > 0, S{t) is a continuous compact operator in . 
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ii) For any Tq £ H'^ , S is a continuous map from [0, t] into . 

iii) {S{ty\t>o possesses an attractor A that is compact, connected and maximal in H^. A attracts 
all bounded subsets of in the norm of , for any r > a — {N/2) — 1. 

iv) If a > {N + 2)/4, A attracts all bounded subsets of periodic functions of in the norm of 
H'', for anyr>a~ {N/2) ~ 1. 

Proof. Items i) and ii) are already proven in lemma 15.61 and lemma 15.51 respectively. 

To verify the rest of items we use results of semigroups and the existence of their attractors (see 
theorem 15. ip . In particular, we need to prove the existence of a bounded subset Bq C -ff*, an open 
subset U of iJ*, such that Bq C U C if^, and Bq is absorbing in U, i.e. for any bounded subset 
B C U, there is a t{B), such that S{t)B C Bq for all t > t{B). This fact is proved in lemma [531 

Item iv) can be checked easily since for a > {N + 2)/4 we have that a > {N/2) + 1 — a and that 
for any Tq G L'^, 

pT 

||AtT||i2 <oo Vr > 0. 



Which together with (|5.5p we can conclude the proof. 



5.2. Weak atractor. We obtain the following estimate for the time derivative of a solution of (|1.5p . 
Proposition 5.8. Let T be the weak solution of (|1.5p obtained in theorem \3. 51 Then, 

^ e LL,(0, (x.; iJ-^) with ^ + 1 = a + ^, (7 e (0, 1). 

Proof. For a smooth ip and < a < 1 we have 



j div{vT)cp <Y^j |A"(i?,t,-'r)||Ai-VI, 



where Rj are the Riesz transforms. Using the inequality (|2.3p with (l/p)+(l/q)=(l/2) and v satisfies 
that ||A°u||Lr < c||AT|jLr, for 1 < r < cx), we obtain 



div(?; • T)lp 



< c||A''r|U.||r|U,||Ai-viu.. 



By interpolation we get 

II' 

Since iJ"/^ c W'P and 
we have 

Therefore, we obtain 



|T||l. < ||T|1^3l|T|l^^ such that s 



N 2 

q 



2N 



P = 



N + 2a-a' 



|A"T|1lp < c|lA2T||i2 and ||T||^_2^ < c|lA2r|li2. 



div(w • T)ip 



<c||r||i^^l|Afr||i+^||Ai-viU2. 



Furthermore, using that weak solutions belong to L°°{0, r; L^)nL^(0, t^H'^) by theorem l3.5| it follows 

— eL/-(0,cx.,i7-i). 

We conclude the proof defining a = 2/(1 + s), cr = 1 — a and using the relations of p, q, a and s. ■ 
By theorem 13. 5l we have that T is a weak solution of (|1.5p that satisfies 

(5.12) T e LL[0, oo; L^) D L^M oo; L') D LL(0, ex.; H^), 

dT N a 

(5.13) — eLLJ0,(X3;ff-^), withy + l = a + -, a e {0,1). 



16 



A. CASTRO, D. CORDOBA, F. GANCEDO AND R. ORIVE 



We observe that Lf^^[0,oo; L'^) is a complete space metric. However, the weak sohitions is not 
closed in this space. We define the space of generalized weak solutions GW{f) formed by the functions 
T G Lf^^[0,oo; L'^) such that are generalized weak solutions of (jl.Sp {T € GW{f )) if T satisfies (|1.5p 
in the sense of distributions and satisfies ()5.12p and (|5.13p . 

Given Ti, T2 G Lp(0, 00; X) {X is a Banach space), we consider the metric 

00 ^ 

(5.14) d(ri, T2) = ^ ^ min (l, |lTi - r2||LP(„,„+i;X)) 

for the set GW{f) (similarly for p — 00). This metric is invariant on Lf^^[0,oo; L'^) (see [2]). Then 
applying classical compactness results (see [22]) for each / e L°°{0, 00, L?), the set GW{f) is a closed 
subset of ifocl^i L^)- Now, we can state the following lemma: 

Lemma 5.9. Given f G L°°(0, cxd, L^) and consider the space GW{f) with the metric defined ()5.14p . 
we have 

i) The set {(T, /)} mt/i T G G'VF(/) and ||/|| < Kq in the norm of the space L°°(0, cx); i^), is closed 
in Ll^[0,oo;L^) x L°°(0, cx); L^). 

ii) The mapping S{t) is a semifl,ow on Lf^^[0,co; L'^) and GW{f) is a positively invariant subset. 

iii) S{t) restricted to GW{f) is compact for t > 0. 

iv) S(t) restricted to GW{f) is point dissipative. 

Proof. 

We verify i) considering a sequence {Tk, fk} e Lf^^i^, 00; L^) x L°°(0, 00; L"^) such that T„ G GW{fn) 
and ||/„|| < Kq. We have that /„ -> / in L°°(0, cx); L^) and T„ -> T in L^^^JO, 00; L^). From the 
estimates of the time derivative of T„ and by the classical compactness results of [23 , we easily obtain 
that T G GW{f) by using the weak formulation of the solutions. Moreover, {(T, /)} is closed. 

Next, we use a smoothing argument and the bound on the time derivative of T as in [21] to show 
the continuity of the semiflow. The positively invariant of GW{f) is immediate. 

The items iii) and iv) are a consequence of the a-priori estimates and the compactness result of 

m- ' 

By the above lemma and using the existence results of global attractors for a point dissipative 
compact semiflows on a complete metric space (see theorem 15. ip . we prove the following attractor 
result: 

Theorem 5.10. Let f E independent oft and a G (0,2). Then, there exists a global attractor 
A, subset of the weak solutions of (jl.Sp . for the semiflow generated by the time-shift on the space of 
generalized weak solutions GW{f). Moreover, A attracts all bounded sets in GW{f). 

We note that the weak attractor is defined in a very weak sense and it gives us less useful information 
than the global attractor in the classic sense. 

Remark 5.11. In the case of time dependent external source f G Lf^^{0, 00; L^) it is possible to extend 
the results of the previous theorem. 

We define ftij) = f{t + t) is the time-shift of f and we consider the hull T-0'{J-) of the positive 
time translates ft with t > of the external force f G T where T C l/^jj^(0, 00; L^) is a set bounded. 
Assuming that T-L'^{!F) is compact with respect to the weak topology of Lf^^{0, 00; L^), then there exists 
a global weak attractor A, subset of the weak solutions of (jl.Sp . for the semiflow S{t){T, f) = {Tt, ft). 
For more details one can refer to the works [2j [21] . 



6. Solutions with infinite energy 

For a divergence free velocity field there exists a stream function V', in the two dimensional case, 
such that 

We shall choose a stream function of the form 

1p{xi,X2,t) X2f{xi,t). 
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Taking the rotational over the equation (|1.2p we obtain 

V X w = d^^vi - da:,V2 = -Alp = -X2dlj = -d^^T. 
Therefore, the function T has the foUowing expression 

r(xi, X2, t) = X2dx^f{xi,t) + g{x2,t) 

where we choose 

1 /■* 

g{x2,t)^-X2 / \\dxJ{T)\\l2,„^^)dT. 



TT 



Substituting the expression in (jl.ip . without diffusion in the two dimensional case, we obtain 
(6.1) dtf. = -dtg - + + g/., 

(here and in the sequel of the section, we denote with subscript the derivatives with respect to x) 
where g satisfies 

t 



1 

TT 

and we define / as 



(6.2) 5(t) = r / \\f-i^)\\h(^.,.)dr 



(6.3) f{x,t)= / Mx',t)dx'. 

J —TT 

Notice that the difference between this system and the one obtained in [J is that this has the property 
of conserving the mean zero value of the initial data. Indeed, integrating equation (|6.ip over the 
interval [— tt, tt) and imposing periodical condition on f^, we have 

dt / Ux',t)dx' = i-Mn, t) + g{t)) / fx{x')dx' . 



Therefore, if fx is a solution of equation (|6.ip and 

/TT 
/,(x',0)dT' = 0, 
-TT 

we obtain 

/■TT 

fx(,x' ,t)dx' ^Q, yt>o. 



6.1. Existence. In this section we prove the following theorem. 
Theorem 6.1. Let ipo e H'^{T) with mean zero value and 

Mo = 'maji(po{x). 

Then, there exist a solution f{x,t) of the equation h6.1\) with initial datum fx{x,0) — (Pq{x) such that 

eC([0,T),i72(T)), 

with T = M^^. 

In order to prove this theorem, first, we add other diffusion term to the equation (|6.ip . Thus, we 
have the following system 

If. 4^ / dtfx = -dtg - ffxx + {fxf +gfx + v{\\fxx\\l2 +g^)fxxx, 

^ ' \ fx{x,0)=M^), 

where g satisfies ()6.2p and > 0. In the next lemma, we prove the global existence of the solutions of 
(El. 



Lemma 6.2. Let ipo e i/^(T) with mean zero value and v > 0. Then, there exist a function f{x,t) 
defined by (|6.3p where fx is solution of the equation |g.^[ ) such that 

/,eC([0,oo),i/3(T)). 
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Proof. We note that if ipo has mean zero value then / has mean zero value. Multiplying the equation 
(|6.4p by fx and integrating over the interval [— tt, tt), we obtain 



L2- 



2 dt 
Therefore, 

Ij^WMlh < +.9ll/.lli^ - H\\f..\\h+9')\\U.\\h- 

Using Gagliardo-Niremberg and Poincare inequalities, we have 

II/.IIl- <C||/.||J.||/,,||J. <C||/,,|U2. 

Hence, 



-Hl/xxlli^-^s'll/xxlli- 



And using the Young's inequality we obtain 
Therefore, 

(6.5) \\M\L^<\\m\\L^-eMC.t), 
and 

(6.6) r \\fxx\\t2dt<C{\\^o\\L-,'y,T), Vr>0. 

Jo 

Taking a derivative over equation (|6.4|) . multiplying by fxx and integrating over the interval [— tt, tt) 
yield 

Hence, 

< < c(||/..||i. +.9ii/..||i.). 

Integrating between and T we obtain 

Jo Jo 



and we can conclude that ll/xajHi^ is bounded for all T < cx). 

Finally, we estimate ||/xa;a;||L2 and ||/2;a:a;a:||L2. Taking two derivatives on the equation (|6.4p . mul- 
tiplying by fxxx and integrating over the interval [—tt, tt] yield 

I d 2 22 2 1 2 2 

2^ll/2:a:a:||L2 + ^^(11/2:2:11^2 + 5 )||/2:a;a;a:||l,2 = - / fxifxxx) + g\\fxxx\\L^ 

< CiT)\\fxxx\\h. 

Applying Gronwall inequality we have that | l/^aa;! |l2 is bounded for all T < oo. We obtain that 
I |l2 is boundcd in a similar form and we conclude the proof. ■ 

In order to prove theorem 16. II we show some estimates, independent of j^, of the global solutions of 
the equation (|6.4p by the lemma [^721 which allows us to obtain the local existence for the equation 
()6.ip . Next, we prove the following lemma. 

Lemma 6.3. Let fx be a global solution of the equation |6'.4| i with initial data ipo and M{t) the 
maximum of fx ■ Then 

(6.7) ^W+^W^T^- 
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Proof. In this proof we use the techniques of article [7] for the control of the maximum of the solution 
of equation (|6.ip . Let denote XM{t) to be the point where reaches the maximum, then 

{M + g)t = M^+gM + i^{\\f,,\\''+g^)f,,{xMit),t) 
< M^+gM<{M + gf. 

Since g{0) = we obtain 
and the proof is finished. ■ 

Proof of theorem l6.ll Multiplying the equation (j6.4p by fx and integrating over the interval [— tt, tt) 
yields 

2di' 

Applying lemma [6?3l and Gronwall inequality we have that | l/xl II^ is bounded for all T < A/(0)^^. In a 
similar way, we can obtain that ||/2;2;||l2 and ||/xa;£c||L2 are bounded for all T < M{0)^^ independently 
of . 

To finish the proof, we consider a sequence of solutions {/'}e>o of the equations 

" dtf^ = -9l - rfL + (Z^)' +57^ + + (5^)')/^ 

/J(x,0) = <^S(x), 

where {'/'olc ^ sequence in H^{T) such that 



rll/.|li2<(M + .g)||/.||i.. 



(6.8) 



/ XX 



lim ip^Q^ ipo e H , 

e— >0 

M'(0) = max/^o(a;) < M(0) = max^o(x), 

and 

||'/5o||//2(T) < ||(po||ff2(T). 

The above estimates provide that 

ll/xllff2(T) is bounded Vr<M(0)~"^ uniformly in e. 
Using the Rellich's theorem we conclude the proof taking the limit e 0. ■ 
Remark 6.4. We shall consider the equation 

(6.9) dtU = ~dtg - lUx + {.fx? + gfx - lyA'^fx, 

which is id.l]) with an extra dissipating term. For this system we have a local existence result similar 
to theorem \ 6.1\ . Moreover, we can construct solutions that blow-up in finite time for a = 1, 2 (see 
below) which show the existence of singularities for DPM with infinite energy. ■ 



6.2. Blovir up. Next, we show that there exist a particular solution of the equation (|6.9p . with v >{) 
and a = 1, 2, which blows up in finite time. 

We consider the following ansatz of (|6.ip and (|6.9p 

(6.10) fx{x,t)=r{t)cos{x), 
then r satisfies 

(6.11) ^^^r{t)[ r^{T)dT + vr. 

dt Jq 

We define the function 

(6.12) Pit) = f r^{T)dT. 



Multiplying the equation (|6.1ip by r{t) we have that (3 satisfies 

^ = 2/3(0« + 2.^. 
dt^ ' dt dt 
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Integrating with respect to the variable t yields 

p'it) - /3'(0) = P\t) - /3H0) + 2v{p - /3(0)). 
Since /3(0) = and /3'(0) = r'^iQ) we obtain 

P'{t) _ . 

r(0)2 + 2vf3 + f3{tY 

If we choose r(0)^ > it follows 



B(t) = Jr(OY - tan (vMOP^^i + aTrta.Ti( — — )) 
Therefore, the function 

(6.13) =r(i)cos(a;), 
is a solution of equation ()6.1|) and ()6.9|) which blows up at time 



- arctanf^^ )) 

" ^r(0)2 -zy2 • 
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